Universal properties of knotted polymer rings 
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By performing Monte Carlo sampling of iV-steps self-avoiding polygons embedded on different 
Bravais lattices we explore the robustness of universality in the entropic, metric and geometrical 
properties of knotted polymer rings. In particular, by simulating polygons with N up to 10 we 
furnish a sharp estimate of the asymptotic values of the knot probability ratios and show their 
independence on the lattice type. This universal feature was previously suggested although with 
different estimates of the asymptotic values. In addition we show that the scaling behavior of the 
mean squared radius of gyration of polygons depends on their knot type only through its correction 
to scaling. Finally, as a measure of the geometrical self-entanglement of the SAPs we consider the 
standard deviation of the writhe distribution and estimate its power-law behavior in the large N 
limit. The estimates of the power exponent do depend neither on the lattice nor on the knot type, 
strongly supporting an extension of the universality property to some features of the geometrical 
entanglement. 

PACS numbers: 36.20.Ey,02.10.Kn 



I. INTRODUCTION 

Scaling hypothesis and the renormalization group the- 
ory led to the understanding of universality in critical 
phenomena, i.e. the fact that systems that look different 
at small scales may share, in proximity of critical points, 
common statistical properties. This concept greatly re- 
duces the variety of possible critical behavior by grouping 
all systems into a small number of universality classes 
characterized by the same critical exponents and criti- 
cal amplitude ratios. In this respect universality justifies 
the study of systems in a given universality class that are 
simple enough to be treated either numerically or ana- 
lytically. 

In statistics of polymer conformations this approach 
has been often followed by modeling polymers as self- 
avoiding walks (SAWs) on discrete lattices such as the 
simple cubic (SC), the face-centered-cubic (FCC) and 
the body-centered cubic lattice (BCC). For these models 
combinatorial arguments and Monte Carlo simulations 
can be efficiently applied to obtain rigorous and numer- 
ical results on the large scale (asymptotic) behavior of 
polymers in good solution Moreover SAWs can be 
mapped to a magnetic system at its critical point and 
studied by renormalization group techniques @, Q • This 
approach has led, for example, to the well established 
results that the number Z(N) and the mean squared ra- 
dius of gyration (Rg(N)) p of A^-steps self-avoiding poly- 
gons (SAPs, i.e. (N — l)-steps SAWs having the two ex- 
tremities one lattice distance apart) behave for large N 
respectively as 



Z(N)^A^N^(l + ^ + ---) 



and 



(1) 
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where the amplitudes A, a, B, b, and the connective con- 
stant \x are non universal quantities that depend on the 
underlying lattice [if, while a and v are universal expo- 
nents related by the hyperscaling equation a = 2 — vd 
with d the dimensionality of the space and v the metric 
exponent. In d = 3, advanced numerical simulations give 
v ~ 0.587597(7) 0, and hence a ~ 0.237209(21), in good 
agreement with field theoretical results [5j . Also the cor- 
rection to scaling exponent A is believed to be universal 
and its best estimate is A = 0.528(12) 0|. 

Notice that SAPs on regular lattices model poly- 
mer rings in good solution with any possible topology, 
i.e. with an arbitrary number of knots tied in. However, 
in most experiments and physical situations such as in 
melts and concentrated solutions, the overall topology 
of a system of ring polymers cannot be changed unless 
the excluded volume interaction is violated. This, for ex- 
ample, has relevant consequences on the temperature of 
the collapse transition of ring polymers, which is found 
experimentally to be different from that of their linear 
counterparts [6J. 

Because of the topological constraint for looped poly- 
mers, the above mentioned mapping to a magnetic sys- 
tem is not valid anymore and, consequently, no field the- 
ory argument can establish the validity of scaling laws 
similar to ([T]) and j[2j . On the basis of previous numerical 
investigations |7Hl0| it is however reasonable to assume 



Z k (N) ~ A k n%N' 
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and 
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where fi k and a k are respectively the connective con- 
stant and the entropic exponent of the subset of SAPs 
with a given knot type k. With this notation k refers 
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either to prime knots or to composite knots given by the 
connect sum of prime knots (roughly speaking, a knot is 
composite if it can be split in two knots located in dif- 
ferent por tions of the chain, which are separated by a 

plane) II1G3. The 

case k — denotes the special case 
of unknotted SAPs (SAPs with trivial topology). 

The scaling relations ([3]) and (0| have been conjectured 
in analogy with ([1} and (J2J) and their validity has been 
confirmed so far only by numerical simulations @, HH 
with the findings /j, k = fi^ < /U, a k = ati]) + m k , Vf. = v and 
B k = Bijj where m k denote the number of prime compo- 
nents in the knot decomposition (Til [l2| . The probability 
of occurrence of a given knot k, 



P k {N) 



Z k {N) 
Z(N) ' 



(5) 



is thus dominated by the exponential decay (^/^,) N at 
large As. In the following we will not deal with this 
aspect but we will focus instead on the scaling beaviour 
of the ratios P k {N) / P^N) = Z k (N)/Z (N). 

It is important to stress that the exponent a k — 
a% + m k is consistent with the recent finding that prime 
knots in SAPs are weakly localized, i.e they occupy on 
average a portion of the ring that scales as N*, with 
t ~ 0.7 .13;]. Indeed, given that each prime knot is weakly 
localized, SAPs with knot type k should, in the large N 
limit, look like as unknotted rings with m k decorated ver- 
tices. These decorations can sit in ~ N positions along 
the TV-step SAP. The partition function of such deco- 
rated chain thus includes a scaling law N nik multiplying 
the partition function of unknotted chains, whose power- 
law part scales ~ N a . Similarly, the average extension 
of each prime component should not contribute, in the 
large N limit, to the overall extension of the knotted SAP 
and one would expect that SAPs with a fixed knot type k 
would share, to leading order, the same metric properties 
of unknotted SAPs Q. 

A more recent numerical calculation of SAPs on the 
cubic lattice has confirmed and improved the above con- 
jectures by suggesting the validity of the following scaling 
laws H 



Z k {N) ~ Z $ {N) x 



N h 



C k 



(0) 



where C k is a coefficient growing with the knot com- 
plexity and factorizable into the elementary contributions 
of each prime component. Roughly speaking C k can be 
interpreted as the elementary entropic cost to tie a given 
knot k in a SAP and its dependence on k can be related 
to the minimal number of steps necessary to build a knot 
type k on the underlying lattice [l4|. Note that, for any 
two prime knots k\ and £/2, the inverse proportionality 
relation A kl /A k2 = C k , 2 / C kl is satisfied. 

One of the aims of this work is to explore more deeply 
the scaling relation ((5]) by extending the numerical inves- 
tigations in [3] to SAPs embedded both on BCC and 
FCC lattices and to look at the dependence of C k on 



these lattices. While the value of C k should depend on 
the three-dimensional lattice considered we expect the 
ratios Z kl (N)/Z k2 (N) to be independent on lattice de- 
tails. This feature has been already suggested in |l5( on 
the basis of a stochastic enumeration of A-steps SAPs 
performed with the GAS algorithm [l|| , an extremely ef- 
ficient method to sample SAPs of moderate lengths and 
with fixed knot type. Note that the conclusions given 
in [ll| rely on an extrapolation to infinity of the data 
obtained for SAPs within the range N < 500, i.e. in a 
region of As where, especially for knotted configurations, 
strong corrections to scaling are expected [17j . 

In this work we extend the range of sampled N up 
to 100000 by means of efficient Monte Carlo samplings. 
We have been able to compute the amplitudes C k for 
different knot types and to confirm, for a wide range of 
N, the universal character of the ratios Z kl (A) j Z k2 (A) 
found in [l5| although with different asymptotic values. 

In addition, by computing the mean squared radius of 
gyration of each SAP, we confirm, within a wide range of 
N, the validity of the scaling law (Q}. Finally, as a mea- 
sure of the geometrical entanglement of knotted SAPs we 
have considered their writhe. Rigorous arguments have 
shown that, for SAPs on the cubic lattice, the mean 
absolute writhe (|WV|) increases at least as rapidly as 
y/~N [l8| and numerical estimates on the same system 
gives (\Wr\) ~ N v with 77 = 0.5035 ± 0.0006 [19]! It is 
then natural to check whether, for SAPs with fixed knot 
type k, the scaling law for the spread of the writhe and 
the exponent r\ are independent either on the knot type 
k or on the underlying lattice. 



II. MODELS AND METHODS 

We perform Monte Carlo sampling of SAPs on the 
SC, BCC and FCC lattices by using the two-point pivot 
moves, a fixed-A algorithm that has been proven to be 
ergodic in the class of all SAPs and shown to be efficient 
in sampling uncorrelated configurations 0, [2(| [2l[ • This 
algorithm allows us to reach values of A up to 10 5 for 
the BCC and FCC lattice. For the SC lattice instead we 
use previous data [3] including also N = 1.5 x 10 5 and 
N = 2 x 10 5 . For the longest chains the samples include 

2 x 10 4 independent configurations. For N — 10 4 this 
number raises to xlO 6 . Typically > 10 5 data are present 
for each A. 

Since the pivot algorithm samples SAPs with any 
topology, it requires a post processing procedure to char- 
acterize the knot type of each configuration. This is can 
be determined, for example, by computing a topological 
invariant such as the Jones or the HOMFLY polynomi- 
als (22j . In general the computational complexity of these 
invariants increases exponentially with the geometrical 
entanglement of the curve and, for very large A, this 
could be the most problematic part of the whole investi- 
gation. To overcome this problem we first smooth each 
configuration by reducing its length while keeping the 



FIG. 1. (Color online) Examples of configurations obtained by the BFACF procedure that reduced their length from N = 100000 
to a few steps, with the hue that follows the monomer sequence: (a) 3i#5i knot on the BCC lattice, (b) 3i#5i knot on the 
FCC lattice, and (c) a very tight configuration with knot type (3i) 4 #4i on the BCC lattice. 



knot type fixed. This is achieved by using a nonequilib- 
rium stochastic scheme based on local moves (a nonequi- 
librium variant of the BFACF method |23l425l |) that trims 
recursively all kinks in the SAP. In addition the algo- 
rithm performs some random local rearrangements to 
reduce the chance that a fast drop in chain length re- 
sults into configurations which it is then very unlikely 
to escape from. This occurs for example when multiple 
prime components are present in the SAP and the re- 
sulting frozen configurations are characterized by having 
these components separated by stretched portions of the 
SAP that become difficult to shrink. We find that in 
most situations the algorithm is able to reach a configu- 
ration close to the minimal N allowed by its knot type 
on that specific lattice. This configuration is then pro- 
jected along a given (approximately) irrational direction 
and the resulting knot diagram is mapped into the cor- 
responding Dowker code [22| that is eventually further 
simplified and factorized into the Dowker codes of the 
prime knots. Finally, we compare each component of the 
original Dowker code against a look-up table of Dowker 
codes of prime knots of up to 11 crossings. In this way 
we have been able to distinguish composite knots with up 
to 5 prime components and with each prime component 
having minimal crossing number up to 11. 

As an example of the efficiency of the simplification 
algorithm we show in Fig. Q] the result of the step reduc- 
tion applied to configurations with N = 100000. The 
first two configurations display a 3i#5i knot (i.e. a com- 
posite knot formed by a prime knot 3i and a prime knot 
5 2 ), (a) on the BCC and (b) on the FCC lattice. One 
can readily see the 3i#5i knot at these reduced lengths: 
the 3i part is on the lower-right side and the 5i part on 
the upper-left one. The configuration in Fig. HJc) shows 
instead a (3i) 4 #4i knot (this a shorthand notation for 
3i#3i#3i#3i#4i). All these examples confirm the effi- 
ciency of the BFACF method in shrinking the chain from 
a very long length to a very short one, even for complex 
composite knots. 



III. RESULTS 
A. Frequencies of prime knots 

The BCC and the FCC lattices have coordination num- 
ber 8 and 12, respectively. It could be expected that more 
neighbors implies more flexibility and thus a more pro- 
nounced tendency to form knots. A manifestation of this 
effect would be that the knot costs Ck were somewhat in- 
versely proportional to the coordination number. It turns 
out that indeed knots are less frequent in SAPs embedded 
on the SC lattice and hence Cf c is larger than its coun- 
terparts for FCC and BCC lattices. However, for the sim- 
plest knots we observe surprisingly that Cf cc < C^ cc . 
According to ^ this is readily seen by determining Ck 
from NZfl/Zk for a simple prime knot. For instance, by 
focusing on in [Fig. Ufa)] an d C l 4 1 in [Fig. [2(b)] it is 
easy to see that these quantities are smaller for SAPs on 
the BCC lattice than on the FCC. If we consider instead 
the ratio Z^/Z^ [Fig. [2(c)] we note that, as N increases, 
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TABLE I. Knot costs C k for the SC, BCC, and FCC lattices, 
determined from averages of the data NZ^/Zk in the range 
5000 < N < 100000 (specific points clearly nonasymptotic for 
the knots with 6 crossings have also been excluded from the 
averages). Errors refer to one standard deviation. The ratios 
^BCCy^sc and zl cc /Zl° in the last two columns seem to 
be independent on the knot k. 
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it approaches an asymptotic value that is independent on 
the lattice. This is true also for frequency ratios between 
other prime knots, as shown in Fig. |3] and in Table |n] 
These results support the previously conjectures hypoth- 
esis that the frequency ratios between prime knots do not 
depend on the lattice in which they are embedded fl5j . 
Our estimates in Table U and in Table [TT| differ however 
from the ones reported in [l5j . A possible explanation is 
that the estimates in [ll| are affected by systematic errors 
due to the small values of N considered whose range fall 
in a region where the corrections to scaling are too strong 
(see the data for N = 400 in Fig. [2fc)) to be neglected. 

Another way of presenting these results is by notic- 
ing that the frequencies of prime knots, whose inverse is 
proportional to the costs C k , have fairly fixed ratios for 
different lattices: according to columns 5 and 6 of Ta- 
ble HI we see indeed that for each knot on the SC there 
are 1.44(3) knots on the BCC, and 1.37(3) knots on the 
FCC (these values correspond to weighted averages of the 
columns assuming a constant knot-independent ratio for 
the two lattices). 



B. Frequencies of composite knots 
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The case of composite knots formed by m k = m topo- 
logically equivalent prime knots k (denoted (k) m here- 
after), was recently discussed in (3] where the entropic 
cost CnAm was determined in a rather simple way in 
terms of the corresponding costs necessary to tie its prime 
components. Since each component is (weakly) local- 
ized, in the large N limit we can think of the chain with 
knot (k) m as an unknotted ring decorated by topologi- 
cally identical objects (prime knots), each placed in one of 
~ N available locations. The entropic cost of each deco- 
ration is Cfc, and being all independent one would expect 
that the full cost C/ k \m is simply given by the product of 
the m elementary costs C k - However, in this picture, be- 
ing the prime knots topologically identical, they can be 
interchanged along the chain keeping the configuration 
undistinguishable. This property gives rise to a combi- 
natorial factor ml that needs to be removed from the 
counting in order to define the relative entropy of the 
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TABLE II. Ratios of partition function of 3i and of other 
simple prime knots k 2 (i.e. Cfc 2 /C l 3 1 ), for the three lattices, 
showing independence on the lattice type. 



FIG. 2. (Color online) (a) Ratio of the probability of unknot 
configurations and of 3i knots, times N x 10 -5 . They con- 
verge to a lattice-dependent constant {Czx x 10~ ) for large N. 
Bands of different colors represent the estimated asymptotic 
values with ± error, see Table [I] (b) The same for knot 4i. 
(c) Ratio of probability of 3i knot and of 4i knot, converging 
to a constant for large N (C^/C^) that does not depend 
on the lattice kind. The horizontal lines mark the interval 
of the estimated asymptotic value, see Table [IT] This panel 
includes also data at N = 400. Note that, while the small 
N deviations from the asymptotic values are due to correc- 
tions to scaling, the fluctuations at larger N's are originated 
by limited statistics. 

decorated chain with respect to the unknotted one. This 
gives 

C (k)m = m! (C k ) m . (7) 

For more general composite knots that include groups of 
different prime knots, the total cost is the product of the 
cost of each group given by ([7]). 

For trefoil knots on the SC lattice is was shown that 
([7]) was fulfilled by an excellent degree of precision [Tij . 
In Fig. 2] we can see that a l so on FCC and BCC lattices 
the data follow the relation Zt^ym /Zq = N m C(k) m with 
C (k)m given by 0. 

The importance of the factorial term m! can be under- 
stood as follows: Since the trefoil is the prime knot with 
the lower cost C k , one would first guess that SAPs with 
the composite knot (3i) m are the most abundant ones for 
TV > Cs 1 . However, when the cost ratio C^/C^ becomes 
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FIG. 3. (Color online) Ratio of probability of 3i knot and of 
(a) 5i knot, (b) 52 knot. The horizontal lines represent the 
estimated asymptotic ratio (Table ITT1 last column) ± error. 



smaller than m, because of the factorial factor m!, config- 
urations with knot (3i) m_1 #4i are more frequent than 
those hosting the knot (3i) m . Since C^/C^ f=s 22 this 
means that the combinatorial entropy loss should start 
to become relevant for a number of prime components 
m > 22, which corresponds to N w (C3.J 22 , i.e. chain 
lengths (for example, N ~ 10 114 on the BCC lattice) 
that are impossible to test numerically. Hence, we can 
safely expect trefoils to dominate the knot statistics for 
realistic polymer rings in good solvent regime. 



C. Metric properties 

By collecting sampled configurations with the same 
knot type we can compute, for example, the mean 
squared radius of gyration (i? 2 (A)) fc of SAPs with fixed 
topology and see if its scaling behavior depends either 
on k or on the lattice in which the SAPs are embedded. 
According to (0} we consider the ratio (R 2 (N)) k /N 2v — 

Bk(l + bk/N A ) and assume A = 1/2 that is close to all 
present estimates [J. In Fig. [5] we report these ratios 
as a function of N^ 1 / 2 for v = 0.587597 [I| and for dif- 
ferent knot types: one can readily see that these ratios 
converge (as N — > oo) to a common value, independent 
on the knot type. This supports the expectation that 
Bk is the same for all knots [7| and confirms indirectly 
the weak localization property of the knotted portion of 
the chain [l]|. The value of Bk does however depend on 
the chosen lattice and is larger on the BCC [« 0.28, see 
Fig.H0>)] than on the FCC [« 0.18, see Fig.^b)]. This is 
somehow expected since Bk is an amplitude. Note that, 
while Bk does not depend on topology, this is not true 
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FIG. 4. (Color online) Log-log plot of ratios of partition func- 
tions Z( 3l )m/Z0 vs. N, for multiple trefoil knots, in (a) the 
BCC lattice and (b) the FCC lattice. The straight lines are 
fits compatible with power laws N m (exponents are in the 
legend). For each lattice, the crossing of these power-laws is 
at a single point of coordinates (C^, 1). 



for the amplitude of the correction to scaling bk whose 
value increases with knot complexity. 



D. Writhe 

The writhe of a curve is a quantity that describe its 
geometrical self-entanglement. A commonly used algo- 
rithm to compute the writhe of a curve goes as follows: 
First, one projects the curve onto an arbitrary plane. In 
general the projection will have crossings that most of 
the time will be transverse, so that, after having estab- 
lished an orientation of the curve, a sign +1 or —1 (de- 
termined by a right hand rule) can be assigned to each 
crossing. The sum of these signs gives the signed cross- 
ing number in this projection. The writhe of the curve 
is obtained by averaging these signed crossing numbers 
over all possible projections. From this definition it is 
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FIG. 5. (Color online) Mean radius of gyration (Rg) divided 
by TV 2 ", with v — 0.587597, i.e. corrections to scaling of (-R 2 ,) 
as a function of TV" 1/2 : (a) BCC lattice and (b) FCC lattice. 
Note that error bars and fluctuations for knots 5i and 52 
at 7Y = 1000 are larger than other ones because of worse 
statistics. 



clear that the main difficulty in computing the writhe of 
a configuration would be the averaging procedure over 
all projections. Fortunately, for polygons on SC, FCC 
and BCC this procedure is enormously simplified by the- 
orems [26l - |28j which reduce the computation of writhe to 
the average of linking numbers of the given curve with a 
finite set of selected push-offs of the curve itself. In our 
calculation of the writhe we made extensive use of these 
results. 

Clearly the writhe of SAPs is, on average, zero and, 
provided we do not distinguish mirror images in chiral 
knots, this is still true for SAPs with fixed knot type. The 
simplest non trivial observable for the writhe distribution 
of SAPs with knot type k is then its standard deviation 
<7fc. In Fig. [5] we report the log-log plot of several <7fc as a 
function of N. Different symbols refer to different knot 
types and it is readily seen that, for sufficiently long N 



FIG. 6. (Color online) Log-log plot of the standard deviation 
of the writhe distribution as a function of /V: (a) BCC lattice 
and (b) FCC lattice. 



all data converge to a common curve that is essentially 
the one for unknotted SAPs. We can then argue that, 
similarly to SAPs in the SC lattice and for all the 
knots considered here <7fc ~ o§ = TV 9 , both for FCC 
and BCC lattices. 

A simple linear fit of the log-log data gives estimates 
% = 0.503(3), D = 0.155(5) for BCC and 7 ?0 = 0.508(2), 
L>0 = 0.1495(14) for the FCC. These estimates are quite 
consistent with the estimate 770 = 0.506(1) found for the 
SC lattice [19j. This is not true, however, for the am- 
plitude Dq, whose estimate in the SC {D$ = 0.1369(7)) 
differs from the ones shown above for the other two lat- 
tices. 

Similar results are found for the absolute value of 
the writhe (exponent 770 = 0.504(4) on the BCC and 
% = 0.508(2) on the FCC), which is expected, given the 
regular shape of the writhe distribution [l9j |. 

Two features are worth noticing: firstly the estimate 
of 7/0 is very close to the lower bound 1/2 proved rigor- 
ously for SAPs with free topology in the SC lattice (18f . 
Secondly the deviation from the scaling law experienced 
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by SAPs with fixed knot type is less pronounced for the 
figure-eight and the unknot (achiral) then for all the oth- 
ers knots that are chiral. This is related to the nonzero 
mean writhe of each of the two images of a chiral knot, 
see [H. 

IV. DISCUSSION 

Since renormalization group arguments are not appli- 
cable to the statistical ensembles of rings with fixed knot 
type, we have resorted to simulations to investigate the 
universal features of the scaling laws for the entropic, 
metric and geometrical properties of knotted SAPs. We 
have confirmed that frequencies of knots in SAPs depend 
on the lattice in which the configurations are embedded, 
On the other hand the ratios between knot frequencies 
are lattice-independent numbers and depend only on the 
knot types involved. Surprisingly, knotted configurations 
are more frequent on the BCC than on the FCC lattice. 
This finding is unexpected because the FCC lattice has 
a higher coordination number than the BCC one, and 
because the knots with the shortest number of steps are 
found in the FCC lattice [HJ. 

We have also supported the conjecture that composite 
knots in SAPs appear with frequencies that, in the large 
N limit, can be inferred from those of the prime knots 
in the knot decomposition. This property was conjec- 
tured for polygons on the cubic lattice [14[ and here is 
confirmed also for SAPS on the BCC and FCC lattices. 
The composite knot cost includes a factorial ml for every 
group of m identical prime knots, which multiplies the 
m-th power of their cost. The factorial term is suitable 



to properly take into account the combinatorial entropic 
loss for identical prime components in the picture where, 
in the large N limit, each prime component behaves as a 
decorating point along the unknotted ring. As a byprod- 
uct one can predict that, for chain lengths N > 10 114 , 
knots with the highest frequency include not only trefoil 
knots but also other prime knots. 

Concerning the metric properties of SAPs with fixed 
topology, our results for all the three lattices considered 
confirm that not only does the mean squared radius of 
gyration share the usual power law N 2v known for the set 
of all SAPs but also the amplitude Bk does not depend on 
the knot type k 0], see (0]). The dependence on the knot 
type is however present at the level of the corrections 
to scaling, which become more pronounced as the knot 
complexity of the SAPs increases. 

Finally we have shown that the large N scaling be- 
havior of the variance of the writhe (or equivalently its 
absolute value) for SAPs with knot type k is independent 
both on k and lattice type and is very close to y/N, which 
is the lower bound proved by rigorous arguments for the 
class of all SAPs on the SC lattice [H]. 

By universality it is reasonable to expect that most of 
the features shown here for SAPs on lattices can be valid 
also for other models of polymer rings in which the ex- 
cluded volume interaction is taken into account. On the 
other hand preliminary studies on the knot probability 
for off-lattice rings have either neglected completely the 
excluded volume interaction (2t| or looked at the knot 
probability amplitudes for few and rather small values of 
N [30( ■ In the near future it would be then interesting 
to extend the analysis reported here to a larger set of 
polymer models. 
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